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ABSTRACT 


The stability of a self -gravitating, non-rotating, plane-parallel, isothermal 
gas layer with equipartition magnetic and cosmic-ray components, immersed in 
a rigid isothermal layer of stars, is considered with respect to waves with 
motions perpendicular to the B e - g e plane, where B, and g e are the equilibrium 
magnetic and gravitational field vectors. The magnetic field and cosmic -ray 
gas hinder gravitational instability, increasing the minimum length necessary 
to produce instability by the factor (1 + a + /3) 1/2 , where a is the ratio of 
magnetic pressure to gas pressure and /3 is the ratio of cosmic-ray pressure 
to gas pressure. 
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I, INTRODUCTION 


In Paper I of this series (Kellman 1972a) we calculated the distribution 
above the galactic plane of an isothermal, plane-parallel layer of gas with 
equipartition magnetic and cosmic-ray components, immersed in an isothermal 
layer of stars. In Paper II (Kellman 1972b) we considered the gravitational 
stability of the gas layer with respect to plane and axially-symmetric perturba- 
tions, neglecting however the effects of the magnetic and cosmic -ray components. 
The result was that the presence of a stellar component increased the radius 
of the marginally unstable state in the symmetry plane, but only slightly. In 
Paper III (Kellman 1972c) the stability analysis was modified by the inclusion 
of a one -dimensional equipartition magnetic field. When disturbances propagate 
across the magnetic field, the minimum length necessary to produce gravitational 
instability is increased by the factor (1 + a) 1/2 , where a is the ratio of magnetic 
pressure to gas pressure. No such simple expression is obtained when disturbances 
propagate along the magnetic field. 

It is our purpose here to include an equipartition cosmic-ray gas in the 
analysis. We consider only the stability of waves (with wave vector k) propagating 
across the magnetic field. Motions are constrained to be perpendicular to the 
B e - g e plane, where B e and g e are the equilibrium magnetic and gravitational 
field vectors. As we stressed in Papers n and HI, the stability analysis may shed 
light (i) on the existence of large-scale structure (1-2 kpc, 10 7 M 0 ) observed 
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in the gaseous component of spiral arms in the Galaxy (McGee and Milton 1964), 
and (ii) on the excitation of density waves in the Galaxy by Jeans' type instabili- 
ties in the gas layer beyond the corotation distance (Lin 1970), ultimately 
leading to the formation of spiral structure. 

U. STABILITY ANALYSIS 

The relation between the vectors B . g » and k and the xyz coordinate 
system is- shown in Figure 1. The basic equations to consider are the continuity, 
momentum, hydromagnetic, Poisson, and heat equation, plus an additional equa- 
tion expressing the fact that the cosmic-ray pressure is a constant of the motion, 
which follows from the fact that the sound speed of the cosmic-ray gas is much 
greater than any other wave velocity considered here (Parker 1966; Field 1970): 


— p + p V’v = 0 
dt * * 

(i) 

V P„ + p *4: v--Lb*VB+-L VB 2 + p V0 = O 

* * dt 477 077 8 

(2) 

^B-Vx ( v x B ) = 0 

(3) 

4ttG (p g + P,) - V 2 (0 g + </>,) = 0 

(4) 

p g =< v L>p e 

(5) 
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( 6 ) 


d 

dt 


Pc-r 


= 0 . 


p, p, 4 >, and § are, respectively, the density, pressure, gravitational potential, 
and magnetic field strength, v is the gas velocity, a first order quantity. 

<v 2 z > is the mean square z turbulent gas velocity dispersion. The subscripts 
g, *, and c-r denote gas, star, and cosmic-ray, respectively. Perturbations 
of the form 


p = p + Ap 

'g 'eg ' g 

(7) 

p = p + Ap = <v 2 > p 

eg ^g tz r g 

(8) 

p = p + Ap 

*c- r r ec - r r c-r 

(9) 

0 g =0 eg+ A0 g 

(10) 

= 4> e + A0 g = 0 eg + + A0 g 

(11) 

B = B +AB 

e 

(12) 


are applied to equations (l)-(6), and terms only to first order in the perturbed 
quantities are retained, with the result that 


— A p + vVp + p V-v-0 

'g 'eg 'eg 


(13) 


— % Ap s + <v t 2 zl > VAp g + P eg ~ v + VAp c _ r + -L V (B e -AB) + p eg VA0 g = 0 


(14) 
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(15) 


d_ 

3t 


AB + v • VB e + B g (V* v) - (B e *V) v = 0 




Ap + v -Vp = 0 

^c-r r ec - r 


(16) 


477G^yO g - V 2 Ac^ = 0. 


(17) 


The subscript 'e' refers to the equilibrium quantities; A refers to the perturbed 
quantities. 

The coefficients of equations (13)-(17) are all independent of t, x, and y, 
enabling us to Fourier analyze in these variables (d/dt - n, d/d x - ikT = 0, 
d/dy - iky), with the result that, written in component form, equations (13)-(17) 
become 


+ ik ^.« v z .♦ (^Tz +Tz p 'i) v - = 0 


(18) 


B 

ik y < V ?x> + n ^eg V y + ik y A Pcr + ik * 4^ ^x + ik y P, t ^g = 0 ( 19 ) 


<vf,> 1 + A „ t ± 4 ,' + n^v, + ^Ap t . r + A ±AB X 


+ — AB — B \p -.0 

47t x dz e eg dz g 


( 20 ) 


ik y B e v y + ( B e-^ + f z B e) v, + nAB x = 0 


( 21 ) 
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(22) 


V 

z 


d 

dz Pec - r 


+ nAp c . r = 0 


-47rGA Pg+ ^ = °- 


(23) 


Equations (19) and (20) are the y and z components of the momentum equation; 
equation (21) is the x component of the hydromagnetic equation. 

As in Paper HI, we assume that 


U c 

O ~ Ip , 

dz eg . e «’ 


(24) 


from which it follows that 


and 



d _ f 
dz Pec * r Pec ' r ’ 


(25) 


(26) 


since B 2 /o and p /p are assumed to be independent of z. Equation (24) 

e 'eg ec-r eg 

is consistent with the presence of a stellar component. Further, one can show 
that 

4-^=-< v L> f (l+“+/3) (27) 

dz e tz 


B e 




(28) 
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and 


1 2a 
p dz y ° g 

eg 



(29) 


where a = B^/8v p eg <v^> is the ratio of magnetic pressure to gas pressure 
and /3 = p ec _ r /p eg < v j z > is the ratio of cosmic-ray pressure to gas pressure. 
We restrict the analysis to the marginally unstable state by setting n = 0. In 
addition, we define the variables e , S, y, and \p by the equations 


e = Ap /p 

8 'eg 

(30) 

S = AB /B 

x' e 

(31) 

y = A Pc . r /<v2 t > 

(32) 

'l> = &t>/<v 2 tz ). 

(33) 


With these various substitutions and restrictions, the system (18)-(23) reduces to 

e + 7/P eg + 2aS + <// = 0 (34) 



One can easily show that 



(37) 
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and defining the variable r by the equation 


t = y/p e g + 2aS, 


( 38 ) 


the system (34)- (36) becomes 


£ + T + 0 = 0 


(39) 


(i|- fa - f/3 ) £+ (i + f ) T+ l ^ =0 


(40) 


- 4 ^i £ + (^-kA^o. 

<v 2 z > \^Z 2 / 


(41) 


We differentiate equation (39) and subtract equation (40), with the result that 


r = (a + / 3) e, 


(42) 


and from equation (39) it follows that 

e = (43) 

(1 +a + / 3 ) 

Equations (41) and (43) are combined to yield a single second order differential 
equation in the variable 0: 


B 2 

B ^ 2 


0 + 


eg 


. 2H 2 (1 + a + (5) P e go 


-k 2 

y 


0 = 0, 


(44) 
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where p ego is the value of p eg at the symmetry plane z = 0. H g is a scale 
height for the gas distribution in the z direction in the limit that a - 0 and 0, 
and has been defined by equation (13) of Paper H: 



< V L> 

8 TrQ P e go ' 


(45) 


The appropriate boundary conditions to impose on equation (44) are 


~yp (z - 0) = 0 


(46) 


limi//= 0. (47) 

| Z | — KX) 

Equation (46) results because y!/ = A0 g / <v* z > is an even function of z; equation 
(47) results because A0 g is constrained to - 0 as j z | — oo . 

Equations (44), (46), and (47) define an eigenvalue problem in the sense that 
only certain discrete values of k y will result in 0's that satisfy (a) equation (44) 
and (b) the boundary conditions imposed by equations (46) and (47). We have 
discussed the nature of this eigenvalue problem in some detail in Paper II. 

Here we need only note that in the limit that a - 0 and /3- 0, equation (44) reduces 
to equation (18) of Paper n. Therefore, to compute the radius r 2 of the marginally 
unstable state (proportional to l/k y ), we simply choose an appropriate value for 
<v? z > 1/2 a*id read off the value of r x from Figure 1 of Paper II, calculated with 
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= 0.025 M 0 /pc 3 , p, D = 0.064 M 0 /pc 3 , and (v 2 ^ 1/2 = 18 km/sec. The 


^ego 

modification introduced by the magnetic and cosmic -ray components is obtained 
by multiplying this r t by the factor (1 + a + /3) 1/2 . The principal result then 
to emerge from this study is that the presence of equipartition magnetic and 
cosmic-ray components tends to stabilize the gas layer against gravitational 
instability, increasing the radius of the marginally unstable state in the 
symmetry plane by the factor (1 + a + j3) 1/2 . In Table 1 we compare the Ledoux 
radius (Ledoux 1951) (p <o = 0, B go = 0, p gc = 0), r x calculated in Paper II 
(p, Q = 0.064 ^ /pc 3 , B eo = 0, p ec _ ro = 0), r x calculated in Paper El (p, 0 = 

0.064 M /pc 3 , B = 3uG, p = 0), and r. calculated here (p, = 0.064 M /pc 3 , 

B eo = 3pG, P ec . ro = 0.50 x 10‘ 12 dynes/cm 2 ), each computed as a function of 
<v 2 z > 1/2 . 

IE. DISCUSSION 

Choosing <v 2 z > 1/2 = 7.5 km/sec and p ego = 1 H atom/cm 3 = 0.025 M 0 /pc 3 , 

l 

values typically quoted in the literature, the Jeans' radius X for an infinite 
uniform gas is 0.638 kpc. All subsequent modifications and improvements to 
the Jeans' analysis tend to increase Kj . Specifically, a self-consistent z distri- 
bution according to the formula p eg (z )/p ego = sech 2 (z/H) increases ^ to ^ = 
0.902 kpc (Ledoux 1951), an increase of 41%. The inclusion of a rigid stellar 
layer with p, Q = 0.064 M 0 /pc 3 and <v^ z > = 18 km/sec increases 

though only slightly, to r = 0.953 kpc (Kellman 1972b), an increase of 6%. 

The inclusion of an equipartition magnetic field increases r t by the factor 
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(1 + a) 1/2 , which for B go = 3/j.G gives 1.10 kpc (Kellman 1972c), an increase of 
15%. Finally, a combined equipartition magnetic \field and cosmic-ray gas in- 
creases r by the factor (1 + a + f3) 1/2 , which for B go = 3/xG and\p ec ro = 

0.50 x 10~ 12 dynes/cm 2 gives 1.31 kpc, an increase of 37%. 
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FIGURE CAPTIONS 


1. The relation between B e , g e , and E and the xyz coordinate system. 
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